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Sign problem of finite-density QCD 
QCD partition function: 


Zaco(T. p) = [DA Det(P(A, pa) + m) exp C Svu(). 
A, —————— 
quark gluon 


Sign problem: Det(D(A, ua) 4- m) Z 0 at uq z 0. 


At T = 0, the sign problem becomes severe at uq > Mz/2 (e.g., Barbour 
et. al. (PRD 56 (1998) 7063)) 
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Baryon Silver Blaze problem: Explain why ng = 0 for ug S my/3 via 


path integral. (Cohen, PRL 91 (2003) 222001) 
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Sign problem of path integrals 


Consider the path integral: 


ja J «eec sta 


e S[x] is real = No sign problem. Monte Carlo works. 
e S|[x] is complex = Sign problem appears! 


If S[x] € C, eom S'[x] = 0 may have no real solutions z(t) € R. 


Idea: Complexify x(t) € C! 
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Lefschetz thimble for Airy integral 
Airy integral is given as 


3 
Ai(a) — EL (5 + 2 


Complexify the integration variable: z = x + iy. 


4 


Integrand on R, and on | 
(a — 1) 
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One-site Fermi Hubbard model 
One-site Hubbard model: 


A = Ufuh, — p(y + fy). 
Fock state gives the number density immediately: 


18 2(eP^ + eBQu-U)) 
bôu D = DF 2e. + PCU)" 


In the zero-temperature limit, 


35 0 08 1 15 4 


(YT, Hidaka, Hayata, NJP 18 (2016) 033002)(cf. Monte Carlo with 1-thimble ansatz: 
Fujii, Kamata, Kikukawa,1509.08176, 1509.09141; Alexandru, Basar, Bedaque,1510.03258.) 
NETA 


Path integral for one-site model 
Effective Lagrangian of the one-site Hubbard model: 
p? 
L= oU + y* [0, — (U/2 + ip + u)] v. 
The path-integral expression is 


B ik d 2 052 
a=a — | ay (14 ef Gem)? 78/20 
27U R eee 


Fermion Det 


Integrand has complex phases causing the sign problem. 
ọ is an auxiliary field for the number density: 
(ñ) = Im(y)/U. 
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Silver Blaze problem for y < —U/2, u < m4/2 


One-site Hubbard model: As BU >> 1 and —U/2 — u > 0, 


U ; 
Det 2, = (» + 3 + ie) E (1 + e &CU/2- n gif)? e» 1. 


The sign problem almost disappears, so that J, ~ IR. 
Finite-density QCD: As 6 > oo and u < m4 /2, 


Det [ya(P4 +m) u] | II (1 + e-805-2)(1 + e- 802) 


= > [. 
Det[u(,-m) nao (+e) +eP) 


(Cohen PRL 91 (2003), Adams, PRD 70 (2004), Nagata et. al. PTEP 2012 ) 
The sign problem disappears by the reweighting method. 
=> Lefschetz thimbles ~ Original integration regions 
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Analysis 


Flows at /U < —0.5 (and u/U > 1/5) 
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Figure: Flow at y/U = —1. Ja 
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Number density: n, = 0 for y/U < —0.5, n, = 2 for u/U > 1.5. 


(YT, Hidaka, Hayata, NJP 18 (2016) 033002) 
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Silver Blaze problem for y > —U/2, y > m,/2 


One-site Hubbard model: At each real config., the magnitude is 
exponentially large: 


Det lo — (» + > + e) =0 (e8 U+#/2)) 


This large contributions must be canceled exactly in order for n = 0. 
Finite density QCD: The situation is almost the same, since 


Det (P(A, uj) + m) _ 
Det(D(A,0) +m) — 


I] ep8 qa), 


Re(AA) «uq 


but ng = 0 for y S my/3. 
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Analysis 


Flows at —0.5 < u/U < 1.5 
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Figure: Flow at y/U = 0 


uU ae 1/2. 


This value is far away from n = Im (z)/U=0, 1, or 2. 


(YT, Hidaka, Hayata, NJP 18 (2016) 033002) 


ru 


Complex saddle points lie on Im(z,,)/U 
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Curious incident of n in —0.5U < u < 0 


We have a big difference bet. the exact result and naive expectation: 


05 O 0.5 1 15 ajo 


This is similar to what happens for QCD and |QCD|. 
p/U = -0.5 & py = Mx /2. 
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Complex classical solutions 


If GU > 1, the classical sol. 
for —0.5 < u/U < 1.5 
are labeled by m € Z: 


u[U x 1/2 


U 
Ea cR By + 27 mT. 


At these solutions, the classical actions become 


SE (B +3)’, 
Re (Sin — Sp) ~ 2m m?, 
Im S. c 2mm (5 4- 1). 


Yuya Tanizaki (RBRC) Lefschetz thimble & Silver Blaze July 27, 2016 € Southampton 12 / 16 


Numerical results 


Results for GU = 30: (1, 3, 5-thimble approx.: Jo, Jo U Ji, and Jo U J41 U Ja ) 
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Necessary number of Lefschetz thimbles ~ 8U/(27). 


(YT, Hidaka, Hayata, NJP 18 (2016) 033002, arXiv:1509.07146[hep-th]) 
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Complex Langevin study of one-site Hubbard model 
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Figure: CL distribution at y/U = 0 


It looks quite similar to Lefschetz thimbles, however the ensemble 
average cannot reproduce their complex phases. 
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TA 
Complex Langevin simulation 
One-site Fermi Hubbard model: 


2.5 


--- exact 
201 HEH w/o phase 
tQ w/ phase 
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(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 


Consequence 
The knowledge of LT helps us to understand failures of CL method. | 
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Summary & Conclusion 


Summary and Conclusion 


e Lefschetz-thimble method gives a suitable framework for 
saddle-point analysis even if S(@) takes complex values. 


e One-site Hubbard model is a nice toy model to play with the 
sign problem. 


e Interference of complex phases among Lefschetz thimbles plays a 
pivotal role for the (baryon) Silver Blaze problem. 
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